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TILT-ROTOR FLUTTER CONTROL IN CRUISE FLIGHT 


Ken-ichi Nasu* 
Ames Research Center 


SUMMARY 

Tilt-rotor flutter control under cruising operation is analyzed. The rotor model con- 
sists of a straight fixed wing, a pylon attached to the wingtip, and a three-bladed rotor. 
The wing is cantilevered to the fuselage and is allowed to bend forward and upward. It 
also has a torsional degree of freedom about the elastic axis. Each rotor blade has two 
bending degrees of freedom. Feedback of wingtip velocity and acceleration to cyclic pitch 
is investigated for flutter control, using strip theory and linearized equations of motion. 
To determine the feedback gain, an eigenvalue analysis is performed. A second, indepen- 
dent, timewise calculation is conducted to evaluate the control law while employing more 
sophisticated aerodynamics. The effectiveness of flutter control by cyclic pitch change was 
confirmed. 


NOMENCLATURE 



fl-t, dj 
Ctjk 




blade aspect ratio 

defined in equations (55) through (61) for k = 1, . . . , 7 
wing aspect ratio 

position of elastic axis from leading edge, nondimensionalized by wing 
chord 

tth, jth generalized coordinate for wing deformation 

kth generalize coordinate for the wing’s jih deformation mode: 

dj , cLj , oy , dy , fl; for k = 1 , . . . , 5, respectively 

position of pitch axis from leading edge, nondimensionalized by wing 
chord 
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C 

Ci 

Qa 


Cu) 

D 

d 

dxjj 

E” 

EIb 

EIc 

EIbw 

EIcw 

Fi 


fm 

f. 

fz 

Gi 


GJw 


G 


Wkj 


h 

I9 

Jb 


matrix, containing virtual work on blade caused by lateral cyclic 
pitch 

matrix, containing virtual work on blade caused by longitudinal cyclic 
pitch 

Theodorsen function 
sectional drag coefficient 
sectional lift slope 

sectional pitching moment coefficient 

position of center of gravity of wing from leading edge, nondimension- 
alized by wing chord 
blade chord 
wing chord 

physical state equation matrix, defined in equations (112) and (114), 
and in appendix B 
drag per unit length 
wing drag per unit length 

physical equation matrix for the nth order time derivative of state 

vector X; defined in equations (112) and (113), and in appendix B 

blade bending stiffness perpendicular to chordline 

blade bending stiffness along chordline 

wing bending stiffness perpendicular to chordline 

wing bending stiffness along chordline 

virtual work on wing caused by blade loading through the tth defor- 
mation mode 

wing torsional moment around the elastic axis (positive nose up) 
wing loading in x-direction per unit length (positive upward) 
wing loading in z-direction per unit length (positive forward) 
virtual work on wing caused by wing loading through the ith defor- 
mation mode 
wing torsional rigidity 

constant virtual work component on the wing, defined by equation 
(105) 

virtual work component on the wing caused by the wing’s jth defor- 
mation mode in term of generalized coordinates Cj and dy for k = 2,3, 
respectively; defined by equations (106) and (107) 
mast height; distance between rotor and the wing elastic Eixis 
unit matrix of order nine 
blade flapping mass moment of inertia 
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thk 



«pp 


'pr 


‘py 


Iw 

^Bkj 


Kh 

Kh 

^Wkj 


Ki 

— cdd^ 
K-u.K — add 

L 

I 

h 

h 

Is 

Iv 

Iw 


constant (independent of generalized coordinates) virtual work com- 
ponent on blade caused by the blade’s j>th deformation mode; defined 
in equation (87) 

virtual work components on blade caused by the blade’s yth defor- 
mation mode in terms of generalized coordinates {ajk); defined in 
equations (93) through (97) for fc = 1, . . . ,5 
virtual work component on blade defined in equation (98) 
virtual work component on blade defined in equation (99) 
virtual work components on blade caused by the wing’s jth defor- 
mation mode in terms of generalized coordinates {qjk); defined in 
equations (88) through (92) for A: = 1, . . . , 5 

pylon pitching mass moment of inertia about the intersection of pylon 
and wing elastic axis 

pylon rolling mass moment of inertia about the intersection of pylon 
and wing elastic axis 

pylon yawing mass moment of inertia about the intersection of pylon 
and wing elastic axis 

wing mass moment of inertia around elastic axis per unit length 
virtual work component on wing caused by the blade’s tth deformation 
mode in terms of generalized coordinates qjk; defined in appendix A 
for A: = 1, , . . , 5 

virtual work component on wing caused by cyclic pitch (cosine term); 
defined in appendix A 

virtual work component on wing caused by cyclic pith (sine term); 
defined in appendix A 

virtual work component on wing caused by the wing’s tth deformation 
mode in terms of generalized coordinates qjk; defined in appendix A 
for A: = 1, . . . , 5 

constant virtual work component defined in appendix A 

constants, associated with first, second order time derivatives of 

wing deformation variables - (uu,,t«ty,pu,), for lateral and longitudinal 

cyclic pitch; equations (116), (117) 

wing semispan 

lift per unit length 

blade lift per unit length 

harmonic component of lift 

steady component of lift 

virtual mass component of lift 

wing lift per unit length 
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M 

Mb 

Mp 

m 

me 

mk 

ms 

my 

Tyi\j) 

N 

Npi 

Nri 

Nyi 

PI 


P 

Pwi 

Pw 

Pdn 



Qni 

R 

Rc 

r 

Tx 

Tz 

S 

Ssk] 



Mach number 
blade mass 
pylon mass 

blade mass per unit length 

wing torsional (pitching) moment around elastic axis per unit length 

harmonic component of wing pitching moment 

steady component of wing pitching moment 

virtual-mass component of wing pitching moment 

wing mass per unit length 

number of blades 

wingtip pitch angle caused by tth deformation mode; equation (44) 
wingtip roll angle caused by ith deformation mode; equation (45) 
wingtip yaw angle caused by tth deformation mode; equation (43) 
virtual work done by the forces on nth blade through the tth defor- 
mation mode 

wingtip deformation matrix, containing Pwi terms, equation (126) 
»th mode of wing torsion deformation Pu, 

wing torsional deformation around the elastic axis (positive nose up) 
nth blade loading in the rotational plane (positive rotational direc- 
tion) 

nth blade loading out of rotational plane (positive forward) 
defined in equations (49) through (53) for fc = 1,. . . ,5 
blade generalized coordinate for jth deformation mode: 

9 j 0 , 4js,Qjc,4jc, for A; = 1, . . . , 5, respectively 

ith generalized coordinate for nth rotor blade deformation 

rotor radius 

radius at airfoil inboard section (rotor cut-out) 
rotor blade radial position from center of rotation 
pylon displacement in x-direction (positive upward) 
pylon displacement in z-direction (positive forward) 
stability matrix, equation (127) 

air loading component caused by blade’s jth deflection mode in terms 
of generalized coordinates (g^); defined in equations (76) through (80) 
for fc = 1,.. .,5 
wing area 

air loading component caused by wing’s jth deflection mode in terms 
of generalized coordinates (cy); defined in equations (81) through (85) 
for A: = 1, ... ,5 


4 



5c 

T 

t 

U 

Uwi 

Us 

y-t 

V 

V 
Vi 
Vr 
V wi 

V^ii 

Vd 

Vin 

Vn 

^out 


Vw 

w 

Wi 

Ww^ 

Wn 

X 

X. 

y 

Zi 

a 

<^G 

Av 

A') 


static mass moment of wing around elastic axis per unit length, Sa = 

TfXxjjCxjD^CLs 

transformation matrix defined in equation (30) 
time 

air velocity at blade section 
fth mode of wing deformation 
induced velocity 
blade tip velocity 
wing upward deformation 

wingtip deformation matrix, containing VV, terms, equation (124) 

airplane cruise velocity 

tth mode of blade deformation Vn 

air velocity at rotor blade section 

tth mode of wing deformation Uu, 

defined in equation (7) 

defined in equation (8) 

downwash caused by wing 

in plane air-velocity component relative to blade caused by wing de- 
formation 

nth blade inplane deformation (positive in rotational direction ) 
out-of-plane air-velocity component relative to blade caused by wing 
deformation 

blade deformation velocity perpendicular to r x n -h K 

wing deformation in x-direction (positive upward) 

wingtip deformation matrix, containing Wwi terms, equation (125) 

»th mode of blade deformation w 

tth mode of wing deformation 

defined in equation (9) 

defined in equation (10) 

nth blade out-of-plane deformation (positive forward ) 

wing deformation in z-direction (positive forward) 

state vector, defined in equation (115) 

wingtip x-displacement caused by tth mode; equation (41) 

wing spanwise position from aircraft symmetry plane 

wingtip z-displacement caused by tth mode; equation (42) 

angle of attack 

geometric angle of attack 

induced velocity 

circulation 


5 



n 

m 

V 2 

e 

6 c 

6n 
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Up 
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correction factor, given by r/irja 
aspect ratio correction factor 

correction factor for spanwise location of vortex trailing 
blade pitch angle 
lateral cyclic pitch 

pitch control input (collective and cyclic) for nth blade 

longitunal cyclic pitch 

wing pitch angle 

collective pitch 

air density 

pylon pitching angle (positive nose up) 

pylon rolling angle (positive clockwise) 

pylon yawing angle (positive counterclockwise) 

inflow angle 

inclination of wing lift 

azimuth angle of nth blade 

rotor rotational velocity 

rotor tth modal natural frequency 

wing tth modal natural frequency 

transposed matrix 

complex variable 


Jl 

a" 


Q Q 2 

time derivative 

Q Q 3 

spanwise derivative 

or radius-wise derivative ^ ^ ^ 
cos{ipn) term; cyclic pitch 
deformation mode i,j 
nth rotor blade 

constant term, i.e., independent of collective pitch 
sin{rpn) term; cyclic pitch 
wing tip 

zero square matrix of order nine 


Superscript 

c 

I 

s 

0 


lateral cyclic pitch 
inverse matrix 
longitudinal cyclic pitch 
collective pitch 
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Subscripts 
B, b 
c 
d 

dd 

i,j,k 

P 

s 

tip 

u 

W, w 
w 

X 

z 


blade 

lateral cyclic pitch 

indicates constant aissociated with first order time derivative 
indicates constant associated with second order time derivative 
indices 

wing torsional deformation 
longitudinal cyclic pitch 
wingtip 

wing upward deformation 
wing 

in K'-constants; wing chordwise deformation 
in i-direction 
in ^-direction 


INTRODUCTION 

The tilt-rotor aircraft is one of the promising transportation systems of the near future. 
It realizes a unique flight capability with its configuration. In the current trend toward 
using active control systems to widen the allowable flight conditions in fixed-wing aircraft 
or in conventional rotorcraft, it is of interest to consider active controls for the tilt-rotor 
aircraft to expand its flight envelope, maneuverabily and versatility. Active control systems 
for rotorcraft have been studied for years, but until recently almost all such studies dealt 
with active control systems to reduce the vehicle’s gust response. It may be said that 
those studies have been aimed at coping with unusual phenomena which occur under 
normal operating conditions. In contrast, the objective of the flutter control studies is to 
enable the aircraft or rotorcraft to fly at higher speeds, to dive more steeply, or to allow 
a reduction of structural weight. In this report, tilt-rotor flutter control using rotor cyclic 
pitch is treated using the governing equations of reference 1. Determination of the control 
law is based on a harmonic method using wing tip deflection as the feedback. The control 
law determined in this way is further examined by timewise calculation using the local 
circulation method, where more precise aerodynamic theories are involved (references 2 
through 5). 
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TILT-ROTOR MODEL AND EQUATIONS OF MOTION 


Figure 1 illustrates the tilt-rotor model considered here, which consists of a can- 
tilevered wing, a pylon, and rotor blades. This model is for airplane cruise mode, which 
is the major concern in this report. Operation in helicopter mode or tranition mode is 
not considered. The wing has one torsional and two bending degrees of freedom about the 
elastic ^lxis. Wing dihedral angle and sweep angle are not considered. The pylon is rigidly 
attached at the wingtip, its direction is assumed to be parallel to the cruising velocity, 
and at the end of the pylon is a hingeless three-bladed rotor. The rotor blade also has 
a straight elastic axis, with elastic lead-lag and flapping degrees of freedom. In both the 
wing and the rotor blade, one of the principal elastic axes for bending coincides with the 
zero-lift- or chord-line. Aircraft motion is not considered. Since the wing is cantilevered 
to the fuselage, wing antisymmetric motion is also neglected (reference 6). 

One coordinate system is assumed fixed to the aircraft (figure 2); consequently, it is 
an inertia system. The X-axis is directed upward. The Y-axis coincides with the wing 
elastic axis. The Z-axis is directed forward. The aircraft cruising velocity V is in the Z 
direction, but the wing-elastic principal cixes do not necessarily coincide with the X- and 
Z- 2 Lxes. The origin of this coordinate system is located at the intersection of wing elastic 
axis and aircraft-symmetry plane. The rotor is located at a distance h forward from the 
wing elastic axis at the wingtip. 

Figure 3 shows the wing and blade deflections. The wing deflection is expressed 
in the X- and Z-direction, and Wu, (positive upward and forward, respectively), and 
torsion about the elastic axis, p«, (positive nose up). The deflected position of the pylon 
root is described as (uu;,<tp, Z., Wu,^tip) or (r*, L, r^). The pylon direction is given by 
yawing, pitching, and rolling angles that are denoted by i/y, i/p, and i/r, respectively, which 
are equal to Pw,tip and — respectively. The rotor rotates clockwise for an 

observer looking forward and the nth rotor blade azimuth angle ipn is measured from the 
vertical position. This reference direction moves with the pylon rolling motion, so the 
blade a.zimuth angle in reference to the inertia system is given by V’n + (figure 4). Blade 
deflection is expressed in terms of in-the-rotational-plane Vn (positive in the rotational 
direction) and out-of-the-rotational-plane (positive forward) bending. 

The blade and wing deflections are described in series of natural modes (reference 7). 

i 


( 1 ) 

( 2 ) 
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Uu, - <li{t)Uwi(y) 

t 

= y^a,(«)lVvv,(y) 

t 

Pw = Yl <ii{t)Pwi{y) 


(3) 

(4) 

(5) 


The equations of motion for the blade and wing generalized coordinates qni and a, are 
taken from reference 1. 

The equations for the deformation of the nth rotor blade are 


Qni + w,-9ni + Vui>r ~ Vot sin 

- ^^Oti'yCOsV'n 


+ Wui>y sin V>n - V^ui>p cos 0n + 2fUVi,i>y cos ipn + 2QWui'p sin tl^n 


= P* 


(6) 

fR 



Vbi = P / mVi dr 


(7) 

Jo 


fR 



V\i= I mrVi dr 


(8) 

Jo 


rR 



W^Ot — Pi rnWi dr 


(9) 

Jo 


.R 



VFii = / mrWi dr 

Jo 


(10) 


where P* is virtual work done by the forces on the nth blade through the tth deformation 
mode . 

For wing deformation the equation is 

"I" Pw i,tip^y) 

+ E E{ 

n=l j ^ 


+ 


+ ^Pwi.up) + 


qnj sin 
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- (^oy-^ + ^ljPwj,tip)^nj COsV>n 


2nVoy( 


+ ^Pw},tip)gnj cos l/)n 


+ 2fi( ^ Voj^Wj,Up)gn 3 sin 
+ ^^(4^oyW"wy,t.p - ^ijPwj,tiphnj cos t(;n 


+ n^ 

= F,+G.- 




9ny sin 


in V'nj 


( 11 ) 


where Fi and G, are the virtual work done by forces on blade and wing, respectively, 
through the tth deformation mode. Natural modes F,-, W(, Uwi, Id^ivt^ and Pwi, and cor- 
responding frequencies a;,- and Wujt satisfy the following equations of motion and boundary 
and normalization conditions. 


Equations of motion for the natural modes of vibration are 
Blade: 


dr^ 


{EIc sin^ 9 + EIb cos^ tf) 


d?Wi 

dr^ 


+ ^ \{EIc - EIb) sin^ cos^ 


-i(i: 


dr^ 

mxdx^^ I — ujfmWi — 0 

dr ' 


*2 


d^V 

{EIc cos^ $ + EIb sin^ $) 


+ I - £/fl) sin^ « COS^ 


dr^ 


i(f 


dW I 

mxdx—r^ ) — m{uj} -I- = 0 

dr ' 


(12) 


( 13 ) 


Wing: 


dj/2 


{EIcw sin^ 0u> + EIbw cos^ 6yj) 


d^Uwi 

dy^ 
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{EIcw ~ EIbw) sin^ dy, cos^ 9^ 


dy^ [ 

^ wiiV^^ ^aPwi} ~ 0 


dy“^ 


dy^ 


{^EIcw cos^ By, + EIbw sin^ ^u,) 


d?Wwi 

dy^ 


+ 


dy^ [ 


j2 tj 

{EIcw ~ EIbw) sin^ By, cos^ 9, 


dy2 J 


- ojlimy,Wwi = 0 


dy 


^GJiv 


- <*>wi{JwPwi + SocUwi) = 0 


The equations for the boundary conditions are 
At the bladetip: 


At the bladeroot: 


At the wingroot: 


Vi" = 0, V/" = 0, Wl' = 0, Wl" = 0 

Vi = o,v' = 0 , Wi = 0 , w; = 0 

Uwi = 0, U{yi = 0, Wwi = 0, = 0, Pwi = 0 


At the wingtip 


Mpfx - — 
ay 

_ ^ 
dy 


{EIcw sin^ By, + EIbw cos^ By,) 
{EIcw - EIbw) sin By, cos B 




Wi 


dy^ J 

d^Wwi] 


dy^ 


+ NMbTx + NMBhPp = 0 


Mpfx- 


^ ■ 
dy 


{EIcw cos^ ^u; + EIbw sin^ By,) 


d^Uwi 

dy^ 


( 21 ) 


_ ^ 

dy 

+ NMet^ = 0 


[EIcw ~~ £/bw) sin cos 


d^Wwi 
' dy^ 


Ipyi^y + {Elcw cos^ d^i + EIbW sin^ tfu;) 




Wi 


+ {Elcw ~ J?/bw) sinful cos 
N 

+ MMsh^i'y + —Ib^v — 0 

£t 


5y2 
d^Wwi 


W r> O 

ay^ 


^pp^p + GJ\\r ~dy — ^ NMsh^i^p + — 0 


/ 0 •> V i 

I pri^r — {Elcw sin Byj -f EIbw cos^ 0^^) 

d^Ww 

- {Elcw - EIbw) sin(?i„ cos^u, — * + NlBi>r - 0 


( 22 ) 

(23) 


(24) 


The equations for the normalization conditions are 
Blade: 


Wing: 


L 


/? 


m(V,2 + ly?) dr = 1 


/ [rriy^ {Uwi + ^Wi) + IwPw + ^EallwiPwi] dy 
Jo 


+ 


Mf {U^i + H'wi) + Ipy (>Vt.)' + IppPii + V (f^Wi)1 , 

Jy=J 

+ nmb \{Uwi + h/v.)" + + W^,] 

L J 

^ [2 {Uwi)^ + = 


(25) 


(26) 


These equations do not have the expressions for the virtual work terms P^, Fi and G^. 
These virtual work terms are directly related to the air loading and can be calculated under 
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various assumptions such as quaisi-steady or unsteady aerodynamics, lifting line theory, or 
lifting surface theory. There are many possible options in calculating aerodynamic loading, 
and some of them restrict the validity of the final equations. For that reason, details of 
the air loading estimation are not given in this chapter and air loadings are assumed to 
be given. Methods to calculate the air loading to complete the equations and to enable 
it to be numerically calculated are given in the following chapters. In the actual aircraft, 
gravitational force also contributes, but this effect is neglected in this report. However, the 
gravitational force depends only on the mass distribution, so it could be easily included. 


Since the blade is rigid in torsion, only lift and drag are considered in the blade 
equation of motions. The direction of lift and drag are determined by the direction of the 
inflow, including the rotor-induced velocity, which also depends upon the method of air 
loading estimation. To avoid additional assumptions, the blade air loading is given by the 
components in the rotational plane (positive in the rotational direction) and out-of- 
rotational plane pzn (positive forward), in the same directions as those used in defining 
the blade deformation. The forces acting in the radial direction are not considered since 
the blade is assumed to neither stretch nor shrink in that direction. Thus the virtual work 
in the equations of motion for the blade is given by 


P'n= f''iPenVi + PznWi)dr 
Jo 


( 27 ) 


As for the case of the blade, wing loading is also given without specifying any aerody- 
namic theory. It is given in cruising direction fz (positive forward) and upward direction 
fx- The wing has a torsional degree of freedom and the torsional moment about the elastic 
axis is given by fm (positive nose up). Then the virtual work G, is 


Gi — 



Uwi + fz^Wi + l^ePwi) dy 


( 28 ) 


The virtual work on the wing by the blade air loading F, is given using the wingtip 
positions Uwi,tip and Wwi,up, and its rotation Pwi,tip- 
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The variable (T) is the transformation matrix from the rotating coordinate system fixed 
to the blade before deformation to the inertia coordinate system fixed to the aircraft and 
is given by 

( 1 -Ur Up \ /cos0„ -sin0n 0\ 

Ur 1 -Uy I I sinV’n cosV»n 0 j (30) 

— Up Uy 1 J y 0 0 1 J 

Blade and wing deformation are expressed by the summation of natural modes for 
each degree of freedom, 


1=1 

(31) 

2 

Y^Viirhniit) 
— 1 

(32) 

1 — 1 
3 

1 — 1 

(33) 

1 — 1 
3 

t= 1 

(34) 

3 

^Pw,(y)a,(<) 
1= 1 

(35) 

3 

1 — 1 

(36) 

1 — 1 
3 

(37) 

II CO 

(38) 

3 

iVpiai(i) 

i=l 

(39) 

3 

^ ^ (0 

(40) 


1 
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(41) 

Z{ = ^Wiytip 

(42) 

Nyi = 

(43) 

Xpi — Pwijtip 

(44) 

Nri = 

(45) 


Even when a simple model is assumed, such as a fully articulated, nonelastic blade, the 
following formulation does not become simple because of the high inflow ratio. Each 
blade is controlled by cyclic and collective pitch; thus 6n and qni are composed of three 
components 


On = 0o{t) + sin Xpn + <?c(0 COS (46) 

Qni = qio{t) + qis(t) sin l^n + gic(t) cos (47) 

Substituting those relations into the equations of motion, the following relations are ob- 
tained 

Blade: 


QtQ+t^i g»o 

-h (gis - 20. qic - qis + w? qis) sin 

-h {qic + 2fl qis - qic + w? qic) cos ipn 

3 

+ + {Qij^j + Qsj^j) sinV'n + {Q\jaj + Qsjaj) cos V’nl 


y=i 

= Pk- (48) 

gij = -F„JV.y+lVoi^ (49) 

Qij = - 

Qij = 2UWuN,j (51) 

= -WuNpi - ^VoiN,i (52) 

Qij = 2muN,i (53) 
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(where Voi, V^it, ^oi, and Wu are given by equations (7) through (10)) 
Wing: 


+ IbNU “ NpiNyj)dj 

;=i 

N 2 f . t 1 

-j- y Ajy 9yo “T Sin [-^iy(vis “ 2n Vic — Via) ^2j 9*0 ^4j 9i0 ~r ^ey VtOj 

n = l j = l 


+ COS xl)n [A\j {(jic + 2fl 9tc) + 9*0 ^5j 9t0 + -^7; 9io] 

+ [Aijiqis - 2H 4ic - Qis) + A\j{qis ~ ^ qic) + qu] 

+ [/Ij,, (s.o + 2fl - n^.c) + A\j (ie + n 9..) + A\, 9i J 

2 

+ ^ s'm 2 tpn[A\j{qic + 21)9,3 - ?ic) + ^3;(9ta ~ 20 9,0 - 0^9,3) 

+ 244 y( 9 ic + 0 9,3) + Asj (9,3 — 0 9ic) + 9 ic + Ajj 9^3] | 


^Fi + Gi 


A\j=V^,Nri+Woj^ 

A),. = -^{Xi + hNpi) + W.jNyi (56) 

A^j = -Voj^Nyi-Wr^N,i (57) 

A‘y = -2n^VojNyi (58) 

4 , =- 2 n^(X, + /iATp.) . ( 59 ) 

A‘^. = 0^ \^{Xi + /liVp.) + VVi>iV„ (60) 

A\^ = ^n^ikVojNy, - RWi.Npi) (61) 


When one compares the terms multiplying the constant (independent of sintpn and 
cos r/j„ terms of equation (48) and the constant term of equation (54), the equations defining 
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Qisi SlXlci. {Zf 2LT6 ol3t;£tlH6d 

3 

qio + <^^qo + J^ Q\j aj = (P* 

y=i 

3 

«i. - 2fJ q„ - (n^ - u,?),., + a, + Qi, ay) = (P‘)^.^ 

y=i 

3 

+ 2 ni, - (n= - wfUi, + a, + dy) = 

y=i 


( 62 ) 

(63) 

(64) 


3 

®t ^wt I B ^^{^yi^pj ^pi^yj)^j 

3 = i 

+ ^ + —^qja + 

y=i 

+ y (^4j + ^^Aj) 9js + y (^^‘sy - 21724*2^) qjc 
+ y (-^ey + ^-^sy ~ ^ -^ 2 y) 9ys 

+ y (-47y — — n^243y) qjc 

= {Fi + G,)q (65) 

where subscripts ( )o, ( )sim ( )cos indicate constant (independent of azimuth angle) 
terms, s'm^n, and costpn terms, respectively, which are considered next. The number of 
blades JV is assumed to be odd and greater than two in these equations, but the calculations 
for an even-bladed rotor can be done in a similar manner. 

The virtual work terms on the right hand side of equations of motion (62) through (65), 
need to be determined next. The terms should be in a linear form to apply the harmonic 
method. The deflection itself, rather than the deflection from the equilibrium state, is 
used in the linearization. Thus, the coefficients of generalized coordinates (i.e., coefficients 
for the first order in Taylor expansion) are different from those at a deflected equilibium 
point. However, the equilibrium value is not correct under the simplifications involved in 
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the air loading calculations, especially using a simple aerodynamic theory. This approach 
was chosen to avoid additional assumptions and to make the following calculations simple. 
Figure 5 gives schematic views of the air loading geometry on the blade and the wing. The 
blade aerodynamic force is calculated using quasi-steady theory. The Prandtl-Glauert rule 
is applied to correct for compressibility effects. The blade lift and drag per unit span, / 
and d, respectively, and the in-plane and out-of-plane aerodynamic forces p^n and pzn are 
written as 


Pen = —i sin <j) — d cos <f> 
Pzn = I cos (f) — dsm<j) 


I = 


1 


2y/l - Af2 


pV^cCiaa 


d — 


Vr = yJV^ + (rn)2 


<f> = tan 

The blade-angle-of-attack equation is 

a = 6 




— (Vin sin (j> - Vout cos <f>) 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 


Vin = Vj + ri'r + {-^x - hOp -t- V Up) sin il^n + {-huy + V Uy) cos (73) 

Vout = IV j + T-z + ruy sin V>n - ri/p cos V’n (74) 

where Vin and Vgut are the air-velocity components relative to the blade caused by wing 
deformation, and is the inflow angle. The induced velocity is neglected, since it is small in 
comparison with cruising velocity, and stability not static deflections is the major concern 
of this study. When Vin and Vout are written in terms of natural modes and Fourier-type 
generalized coordinates, the following relations are obtained 


Vins'm<i) - Vout cosd> 

2 3 

= ^ Setjqjo + Swijdj 



( 75 ) 


18 


Ssij = Vj sin<^ — Wj cos<j) (76) 

SB2j = Sbij (77) 

Sb3] = (78) 

Sb 4} = Ssij (79) 

Sb 5J = ^Sbij (80) 

<5^1; = sin (j> — Zj cos <f> (81) 

Sw 2 i - - {Xj + hNpj) sin <f> - rNyj cos <f> (82) 

S\v3j = ^ ^pj sin (f) (83) 

S\V4j = —hNyj sin<^ + rNpj cos <j) (84) 

S\v5j — V NyjS\n<f> (®^) 


The relations (66) through (85) give blade and wing air loading in terms of the gen- 
eralized coordinates. Substituting equations (66) and (67) into (27), the virtual work on 
the blade caused by the blade air loading is given as follows: 


Pn -^BO + ^ ^Blj Qjo + ^ Iwij 


j=i 

2 


J = 1 


+ I ^ ^BSj^jc ^W2j^j ^WSJ^j + Is^s 1 sin tl^n 

\j=l j-1 j=l j=l J 

( 2 2 3 3 \ 

^B4j^jc + ^ ^W4j^j ~ ^ j COS Ipn (86) 

j=l j=l j=l j=l J 


I'bo = r - IV.cos,^)* 

1 

--P cCdVl {Vi cos <I) + Wi sin <t>) dr (87) 

2 Jo 

1 1 

^Bij = -nP^la / VRcSBij{Vism<i>-WiCos<t))dr (88) 

^ Jo \1 ~ 

1 1 

lB 2 j = -7.P C'/a / -^=^=VRcSB 2 j {Vi sin (f) - Wi cos 0) dr (89) 

^ Jo V 1 ~ 

^3, = -\pCia r ~^J=VRcSBz,{Vi3m<i>- WiCOs<t,)dr (90) 

* yo V 1 ~ 


19 


1 1 

lB4j = -^=^VRcSB4j{V,sin<t>-WiC0S<f>)dr 

4$; = -j===VrcS Bsj {Vi s\n <f>-Wi COS <l>) dr 

1 1 

4iy = -^^==V«c5vviy(V^.sin<^-VyiC0s<^)rfr 

1 1 

= - 2 ^ v/l - ~ 

1 /*^ 1 

7^3^. = --pCla J 

1 /* ^ 1 

^W4y = "2^ J^ v/l- ~ 

^W5y = - ^ yi - ^’’ 

Ic = -7;P^ic,f ]:V^{Vi sin (f>-Wi cos <f>) dr 

I Jo ^ 


(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 

The virtual work on the wing done by the rotor-blade air loading, given by equation (29), 
is 


7*c = Pn 


N 

=£/. 


Uwi,tip 
0 

y^Wi,tip , 
W(vi,Up 


0 

(T) I pen 
Pzn 


+ Pwi.tip 1 {T) 


-f/' 

^Wi,tip 


r 

Vn 
h-\- Wn 


0 

X 1 pen 
. pzn 


dr 


^ O 

= XI -^vviy^j + X ■^vv2j°y 


y=i 


y=i 


+ X ^Blj9j0 + X 

y=i y=i 

2 2 

+ X ^Biy9y0 + X ^B2j^j‘ 


y=i 


y-i 
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2 2 

+ ^B4jQjO 

j = \ j=l 

2 2 

+ ^ + Y1 ^B6j9jc 

j=l j=l 

+ + Kl, 

Details of Kg^ given in appendix A. 


(100) 


The wing aerodynamics are more complicated because of the pitching motion. When 
some aerodynamics theories are referred to as quasi-steady, the reference often means 
the usage of wing velocity relative to air, neglecting the effect of shed vortices. When 
pitching motion is introduced, wing velocity relative to the air varies along the chord. 
Aerodynamic force becomes dependent upon the point used to represent the wing motion. 
It is not appropriate to take the elastic axis as a representative point, since the air forces 
are not directly related to blade elasticity. In steady or unsteady aerodynamics (such as 
lifting-surface theory or Theodorsen’s theory) the three-quarter chord line is an important 
point to represent the boundary condition. In this formulation the boundary condition 
for the angle of attack is at the three-quarter chord line, and the effect of shed vortices is 
disregarded. The wing aerodynamic force is given by lift 1^, drag and pitching moment 
around the elastic axis mg 


I ^.t/ 2. Ffl I _ 

Iw — 2^^ X \yw Pm 

* ^ ^ T/i . + Wiy 

^ “t" Pm 

TYle “ 7\P^ 0.25)c^£> 

It 


(0,75 ^e)^wPw 
V + tbxju 
( 0.75 
V 


1 


+ 


y2 \ 


( 101 ) 

( 102 ) 


where Cm is the pitching-moment coefficient about the aerodynamic center (which is as- 
sumed to fall on the one-quarter chord line), and CgCu, is the distance between the leading 
edge and the wing elastic eixis. Downwash caused by the wing itself is denoted by Vd and 
is given using the geometrical angle of attack 6^ and the wing-aspect ratio Aw 


Vd 


“2- Aw + \/4 + 

2+y/rr^ 


( 103 ) 


The inclination of lift is expressed as 0^, = VdjV . When this inclination of lift is assumed 
unchanged to be, the aerodynamic force on the wing is given in the X and Y directions 
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as {Ucos<f>uj~0.5pV‘^Cu;CdSin<j)u^, -lwsm<i>u, +0.5pV'^CwCd<^os<f>u;). The virtual work on 
the wing done by the wing aerodynamic forces is given by 


• L / lujCos4>u, - 0-5pV^CujCdsin<t>t 


rL h XD w oaxi v^ty 

G, = / —l^s\n<i>w + Q-^pV^c^CdCos<f>w 1 | ^Wi j dy 

Jo i 


Uwi 

Wwi 

Pwi 


= *^VV1 + ^ ^W2j ^ ^\V3j 

j=l j=l 

G\y \ — f Cu,G(q.(V $uJ ~ Vd)‘ 

Jo ^ 

\Uwi cos - Wwi sin ^xu + (^e - 0.25)cu,Pwt] dy 

/•L ^ 

- / -/)V^Cu,Gd(f/wtSin^u> + W^wicos0u,)dy 

Jo 2 

+ IpVhlCmPwidr 

Jo ^ 

f^l 

^W 2 j ~ J Pw j* 

\Uwi cos <i>xxf - Wwi sin (f>xjt; + (og - 0.25)cyjPwi] dy 
^W3j ~ [ 'ZP^w^la[~^^W j “t" (0*75 — ^wPw j ~ ^d^Wy]* 

Jo 2 

[(7vvtCOs</>u; — W^vvtsin</> u> (ce — 0.25)ciyPvvt] dr 


(104) 


(105) 


(106) 


(107) 


After the virtual work is written in terms of generalized coordinates, and when one 
compares the terms multiplying the constant (independent of azimuth angle V’n)? sin^n 
and cosxpn terms, the following equations are obtained for the blade equation of motion 

3 2 3 

Qio + 9to + ^ Qij dj - ^ pBijQjo - ^wijdj = Pbo (1^®) 

y=i y=i ;=i 

3 3 

Qis + (t<^t* ~ ^^)Qis - 2fl Pic + ^ Q^jdj + ^ Qzjdj 

j=i j=i 
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2 2 



j=i i=i 

3 3 


J=1 J=1 


3 3 

9tc + (w? - f^^)9tc + 2Q g,a + Q\j^j + ^ Qsj®; 

y=i j=i 

2 2 
~ ^BAj^jc - y^^BSjQja 

y=i y=i 

3 3 

“ X!/ ^WAj^j ~ Y1 ^WSjO'j = Ic^c 

j=\ y=i 


( 109 ) 


(110) 


In the same way, by substituting equations (100) and (104) into equation (54), the 
following equation is obtained from the constant term of the wing equation of motion: 

3 

+ ^wi^i ~ ^pi^yj) ~ ^W3j ~ ^Wlj] 

J = 1 
3 

“ XI(^VV2j + ^W2j)°‘j 

i=i 

2 2 2 
+ X^ — X^ ^BljQjO + X^ ^BAj^jO 

j = l j = l y=i 

^ N ^ f N \ 

■*■ X^ ~^-^2j^ia + XZ ( ~ ^B2j ] 9js 

J=1 J=1 ^ ^ 

2 

X^ ~ ~ 9ja 

+ X^ y ^3j9yc + X^ ( y ^5; “ N?lA\j - j Qjc 
+ X^ ( ~ ~ “ ^B6j^ Qjc 
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= K*g6e + + Kq (111) 

Final equations of motion for the generalized coordinates in equations (108) through (110) 
are written in the following form: 


E^X + E^X E°X 

1! 

0 

o 

s <?c)* 

(112) 

E" = 

/Ef,. 

jpn 

^2,1 

xpn 

^1,2 

jpn 

^2,2 

• • • Ei g ^ 
. . . E2^g 


(113) 



j7»n 

^9,2 

• • • Eg g j 



D = 

<£>1,1 

£>2,1 

£>1,2 

£>2,2 

£>1,3 \ 
£>2,3 



(114) 


£>9,1 

£>9,2 

£>9,3 ^ 




O 

II 

9lC? 

920, 92s, 92c, «1, 02, O3)* 

(115) 


Details of the elements and Ay are given in appendix B. 

DESIGN OF CONTROL SYSTEM 

One approach to suppress coupled wing/rotor self-exited oscillation is active pitch 
control of the rotor blade (references 8-11). The tilt-rotor aircraft has other control sur- 
faces, such as a flaperon. However, only the pitch control using a conventional swashplate 
is considered here. The control system is a simple feedback system which uses wing- 
deflection feedback. The possibility of flutter control is investigated using a harmonic 
balance method. For simplicity, air loading is calculated here using a two-dimensional 
strip theory for both the wing and the rotor blade. In the following analysis, the gravita- 
tional force is neglected. When blade pitch control is introduced, the natural frequencies 
and corresponding mode shape of the system vary as a function of time. This change is 
not considered and the natural frequency and mode shape are assumed to be valid for a 
specific flight condition. 

This formulation uses generalized coordinates. What actually can be measured in 
flight is not a generalized coordinate but the actual deflection or its derivative, which is 
given by an infinite series of generalized coordinates. There are two choices in introducing a 
feedback control system: one uses a generalized coordinate as input for the control system, 
and the other uses an actual deflection or a rate of deflection. The latter is chosen here. 
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The wing deflection speed and the acceleration at the wingtip are used for the feedback 
input. Cyclic pitch control is given by 

0 s = Kxisd^Wjtip f^usdd^Wjtip 
+ Kxifsd^Wftip ^wadd^Wytip 
J^psdPwytip “1“ J^psddPwytip 

0c ~ f^ucd'0'Wytip J^ucdd^Wytip 
"I" f^wcd^Wytip “I" f^wcdd^Wytip 
“H KpcdPwytip ”1“ ^pcddPwytip 

Using this control law, the final equation is modified to 

E2 <=X + E^"X + E°‘=X = D(^o 0 0)* (118) 


E^-^ = E^ - + K^cddB’^)y 

- {Kr,sddB^ + K^cddB<=) W 

-{Kj,sddB^ + Kj>cddB^)P 

E'" = E' - {K^sdB^ + KvcdB*^) V 

- (K«,.dB^ + K^cdB^^) W 

-(Kp,dB* + /fpcdB‘=)P 

j,0c _ jjO 


(II 

Is 

Is \ 

0 

0 

0 

0 

0 

0 

r2 


^s 

0 

0 

0 

0 

0 

0 

K\ 

Kh 

Kh 

Kl 

Kl 

Kl 

\Kl 

K% 

K%) 


( 119 ) 


( 120 ) 

(121) 


(122) 
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V = 


w = 







(Ih 

Ih 







0 

0 

0 







0 

0 

0 







Ih 

Ih 

Ih 




= 

= 

0 

0 

0 







0 

0 

0 







Kh 

Kh 

Kh 







Kh 

Kh 

Kh 







\K% 

Kh 

Kl) 


fo 

0 

0 

0 

0 

0 

Uwi,tip 

ywi,tip 

0 

0 

0 

0 

0 

0 

Uw2,tip 

V\V2,tip 

Vo 

0 

0 

0 

0 

0 

yw3,tip 

VvV3,ttp 

/O 

0 

0 

0 

0 

0 

Wwi,fip 

Wwi,tip 

0 

0 

0 

0 

0 

0 

Ww2,<tp 

Ww2,ttp 

VO 

0 

0 

0 

0 

0 

Wws.ttp 

Ww3,ttp 

/O 

0 

0 

0 

0 

0 

Pwi,tip 

Pw 

l,tip 

0 

0 

0 

0 

0 

0 

Pw 2 ,tip 

Pw 

2, tip 

Vo 

0 

0 

0 

0 

0 

P\V3,tip 

Pw3,tip 


VWl, tip 

Vw2,tip 

Vws,tip 

^Wl,tip 

W^W2,tip 

M^W3,tip 

Pwi,t 
Pw2,tip 

PwZytip 


(123) 


(124) 


(125) 


(126) 


The stability of the system, whose motion is given by equation (118), is defined by 
the eigenvalues of the following matrix S. 

(127) 

where I 9 , and O 9 indicate the inverse matrix of a unit matrix of order nine, and 
a zero square matrix of order nine, respectively. 



TIMEWISE CALCULATION 

To make a transient numerical calculation for the tilt-rotor’s coupled-wing and blade 
deflections, and to evaluate the control law given in the previous section, a method is 
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given for independent timewise calculations. Determination of the control law is based 
on a rather simple theory. A more sophisticated aerodynamics theory and two or more 
natural oscillatory modes for each degree of freedom are employed in the transient analysis. 

Most of the aeroelastic studies conducted previously have been using quasi-steady and 
strip theory. In this study the aerodynamic force is calculated on the basis of an unsteady 
theory. The calculation is based on Weissinger’s concept through a vortex lattice method 
for the wing and on a local circulation method for the rotor blade (both with modifications 
to take the shed- vortex effect into account). Unsteady aerodynamics is a very complex 
problem even for a fixed wing with infinite aspect ratio. To get a precise solution requires 
large amounts of computer time. As for this study an efficient timewise calculation of 
tilt-rotor motion was desired, the following formulation aims at a rapid calculation rather 
than an exact solution. Arbitrary wing or rotor-blade motion is treated as a superposition 
of harmonic motions for a short interval to make the calcination simple. Then, for each 
harmonic motion, wing theory in harmonic motion is applied. 

Wing lift lu) and pitching moment me consist of three components: steady components 
Is and ms, harmonic components Ik and mk, and virtual-mass components /i, and m^: 

Iw = h + ^ h + Iv ( 128 ) 

k 

me = ^ mfc + mt,. (129) 

k 

As the shed vortex distribution depends upon the history of the wing motion the wing 
deformation and its derivatives are expressed as a series of harmonic components, using 
their previous values: 


k 

(130) 

Xl\D ~ “ 1 ” ^ ^ “^wk^ 

k 

(131) 

VJ\jD ^ "^wk^ 

k 

(132) 

k 

(133) 

Pw ~ PwO ^ ^ Pwk^ ^ 
k 

(134) 

Pw ~ PwO ^ ^ Pwk^ 

(135) 


k 
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where variables with a bar represent complex numbers. 

The boundary condition at the position bc^ from the leading edge is 

ti = -titu + Pw{b - auj)cyu + (V + w W )(^w + Pw) (136) 

When the deflection and its time derivative are assumed to be small, the equation is 
linearized as 

U — ^wj^w ^ “ 1 “ Pw") ” 1 " 

The steady lift component is calculated through Weissinger’s concept. At the three- 
quarter chord line, the downwash velocity is given by 

Uo,75 — ^ tv ”1" Pw) ”i~ (f^^) 

Since this value is dependent upon time, it is not directly used in the boundary condition 
for the air loading calculation. The steady component of the air loading is calculated using 
the constant term of this air velocity as the boundary condition in a vortex lattice method 

tic — tiujQ ”f“ Py,o(0.75 ^tv'^^w Pwo) ti^wO^w (f^^) 

Steady pitching moment around the elastic axis rris is attributed to wing airfoil shape and 
is given using the pitching-moment coefficient for a steady airfoil 

~ /^(0.25 — o,i,)cu; 

The steady lift given in this way is accurate in the sense that the shed vortex effect does 
not exist and that the trailing vortex effect is included. 

The unsteady components (harmonic components and virtual mass components) are 
more difficult to estimate. Even under the assumptions just given (sinusoidally oscillating 
air loading), there is not an easy and versatile calculation method which properly includes 
shed vortex and trailing vortex effects. As stated previously, this study aims at the efficient 
timewise calculation of tilt-rotor motion. Air loading /fe, caused by the A:th harmonic 
motion is calculated using a two-dimensional theory. Aspect ratio and spanwise location 
are taken into account by multiplying /jt by a factor q: 

Ik = rfCiapV ~C (-Uwk + Vpwk + (0.75 - aw)ctuPu,jt + ^wk^w) (141) 

fnk = -lk(0.25 - au,)c«, (142) 
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where C (u>a;Cui/ 2V') is the Theodorsen function, and the factor Tf is given by 


1\ = V 1 V 2 

(143) 

2 + \/4 -h AIj 

(144) 

4 / / w \ 2 


(145) 

II 

3 

(146) 


The variable r/i is the correction for finite aspect ratio Axv, and r }2 is calculated under an 
assumption that this air loading is distributed elliptically along the span. The air loading 
caused by the virtual mass motion is also given by two-dimensional theory with the same 
correction: 


Iv — + ^ 3 ) 

(147) 

rriv = r){m 2 + m 3 -h m. 4 ) 

(148) 

h = />7t (^y) + (0.5 - au,)cu,pu.) 

(149) 

I 3 = pT^ (y) y?w 

(150) 

m 2 = -h (0.25 “ Cw 

(151) 

7723 — /3(0.T5 ^w')^w 

(152) 

1 / ^ ^ 

(153) 


Drag is calculated using a drag coefficient 

d = -pV^CyuCd (154) 


These equations give the lift, drag, and pitching moment, but not the directions of 
lift and drag. In the Ccise of lifting-surface theory, the most accurate expression of lift 
and induced drag is the sum of integration of lift along the chord line and leading-edge 
suction. Because of the simplification, where the chordwise distribution of bound vortex is 
not taken into account, this expression is not available in the present method. The most 
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suitable form for this simple case is that the lift is inclined backward by an amount, equal 
to the angle between the free flow and the lifting line (the same expression as that in simple 
lifting-line theory). It is sometimes discussed in the term of induced drag, and it is known 
that the estimation of induced drag is difficult even in a complete lifting surface theory. 
The lift and drag values are converted into air-load components in x and z direction {lx 
and Iz) 

lx = lcos(f>u) — ds\n<)>^ (155) 

Iz = —lsin(f>^JJ - dcos<p^ (156) 

The calculation of lift-inclination angle <f>w is given in appendix C with the formulation of 
a vortex lattice method. 

The blade air loading calculation is similar to that of the wing, except that it is based 
on the local circulation method instead of a vortex lattice method and that it lacks the 
torsional degree of freedom. The details of the local circulation method are described in 
refernces 2 and 3. Here the method used in the unsteady modification is stated. Among 
two dimensional theories of harmonic air loading, Theodorsen’s and Sears’ theories are 
widely used. Blade deflection in two directions, in plane and out of plane, is treated just 
like the wing case. Apart from the deflection, the blade encounters the induced velocity 
field made by the returning vortex; hence the local circulation method is better explained 
by the Sears’ theory. However, as stated previously, because this model does not have 
a pitching degree of freedom, and because the Sears function can be approximated by 
the Theodorsen function when reduced frequency is small, this remaining induced velocity 
is included in the boundary condition on the blade element. Harmonic air loading is 
calculated using Theodorsen’s theory in this study (reference 12). The air loading caused 
by virtual-mass effects is added afterward, because the relation of induced velocity and 
circulation is given by the quasi-steady assumption; thus, it does not contain virtual mass 
or acceleration terms at all. In this calculation, variation of pylon direction is assumed to 
be small, so rotor-blade air loading is calculated using the theory for axial flight. 

The components of blade lift are steady lift, harmonic lift, and virtual mass lift 

lb — Is + ^ 

k 

The blade motion relative to the air is composed of cruising velocity, rotating velocity, 
induced velocity, and blade- and wing-deformation velocity (figure 6). These velocities 
give the geometrical angle of attack at the elastic axis as 

l)p + Us COS (f) 

OLG = 0 - 4> 
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. V 

<p = tan — — 

rn 


U = \/(rn)2 + y2 


(159) 

(160) 


When Vp is written in terms of in-plane and out-of-plane deformations equations (158) 
becomes 


a A 

ac = 0 - 4>- — 


(161) 


Ut = {Ua — w — Wfju — rUy sin Ipn + cos Ipn) cos <i) 

H- [v -f- ri/f + {—iiw — hvp -t- Fi/p) sin V>r» + {~huy -f Vi/y) cos ipn] sin <l> (162) 

This geometrical angle of attack is written in Fourier series, based upon its past values 

0!G = aco + ^ (163) 

k 

For the steady component (Xao, the conventional local-circulation method with the Prandtl- 
Glauert rule is applied, and the following equation is satisfied at each control point: 




2y/i 


(«<,„_ i^A„) 


= pV^Al 


(164) 


where At; and A 7 are the induced velocity perpendicular to local airflow and the circulation 
of the hypothetical wing, respectively. 

The sinusoidal component of the air loading is given by two-dimensional theory, with- 
out considering the effect of compressibility, but with corrections, similar to the wing-air 
loading case, for aspect ratio and spanwise location: 


fife = ^pU cCiaC ^ OCGke * 

(165) 

T) = T}iT]2 

(166) 

2 + y/4 + A^ 

(167) 

V2 = -V^l - 

7T 

(168) 

R — 

A= ^ 

c 

(169) 

^ 2 r — (i 2 + Rc) 
^ - R-R^ 

(170) 
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Theodorsen’s theory treats both the pitching motion and the heaving motion of a wing. 
The relative velocity of the blade to the air varies along the chordwise position because of 
pitching motion. Since the rotor blade model does not have a pitching degree of freedom, 
virtual-mass terms are assumed to be related to the acceleration of the blade elastic axis 
neglecting blade pitch motion. 


= pn {w cos <l> — V s'm 4>) 


(171) 


This equation concludes the air loading calculation. The virtual work required to 
complete the equations of motion can now be calculated by equations (27) through (29). 
Timewise integration is carried out using a Runge-Kutta-Gill method. By conducting 
timewise calculations both with and without the feedback system, the effectiveness of the 
feedback system can be evaluated. 


NUMERICAL RESULTS 

The methods of control law determination and the timewise calculation just given 
are applied in this section. The specifications of the model aircraft used in numerical 
calculations are given in table 1 and in figures 7 and 8. Blade pitch angle is chosen to 
give the thrust of 7,OOON at the design advance ratio V jQR = 0.7 and such that the 
thrust changes proportionally to the square of cruising speed. Trim thrust and pitch angle 
determined in this way are given in figure 9. 

Figure 10 presents the change of the real part of the eigenvalues as a function of aircraft 
cruise velocity. This figure shows that this model remains stable for advance ratios up to 
1 . 8 . 


To investigate the feedback-control efficiency, the stiffness of the considered tilt-rotor 
model is reduced. Because flutter peculiar to the tilt-rotor aircraft is the major concern in 
this report, the wing stiffness was reduced to one-eighth of the baseline model for all the 
degrees of freedom in the wing, while the rotor-blade stiffness is assumed to be unchanged. 
Natural frequencies and mode shapes employed in the eigenvalue analysis are given in 
Figures 11 through 15 for the new rotor/wing model. These modes are attributed to 
blade-in-plane, blade out-of-plane, wing-upward, wing-forward, and wing-torsional degrees 
of freedom, respectively. Figure 16 gives the real part of the eigenvalue as a function of 
advance ratio V /QR for the reduced stiffness model. An eigenvalue will become unstable 
when the advance ratio exceeds about 1.5. 
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For further numerical results, an almost critical case in which the advanced ratio 
is equal to 1.5 is chosen and the design of a control law is investigated. In the previous 
section, input to the feedback system as a function of wingtip upward and forward bending 
and torsion was described. Using all of these parameters is not necessary in the present 
case. When selecting one or more of the parameters, the eigenvector and mode shape 
corresponding to the unstable mode is of greatest importance. The two eigenvectors which 
have the greatest real part at the advance ratio of 1.5 are given in figure 17. Among the 
elements of the eigenvector, only the elements which have absolute values of more than 5% 
of the largest element value are shown in this figure. The critical element is attributed to 
the wing’s second lowest natural frequency and the other element to the the wing’s third 
lowest frequency. Wing second and third modes correspond, in turn, to forward deflection 
and torsion, respectively, as seen in Figures 14 and 15. To stabilize or to augment the 
stability of these modes, feedback of dominant deformations are expected to be effective. 
As the simplest case, paying attention to the least damped mode only, velocity of the wing 
forward bending is chosen as an input for the feedback system, which means all the 
feedback gains except for Kwad and K^ocd are equal to zero. 

Figure 18 shows the effect of feedback on the real part of the systems eigenvalues. Only 
those eigenvalues that have a real part greater than —1.0 are shown among the eigenvalues, 
since highly damped modes need not be considered. The imaginary parts (frequency) of 
the eigenvalues change at most 2% when the feedback gain changes from 0.0 to -0.05. 
The eigenvalues without feedback are given along each line. The dominant, generalized 
coordinate corresponding to each eigenvector is also given along each line in parentheses. 
All the modes are simultaneously stabilized when feeding back the cosine term of cyclic 
pitch. Some modes are stabilized yet others are destabilized in case of feedback with 
the sine term of cyclic pitch. Feedback of the sine term of cyclic pitch works about as 
well as that of the cosine term as far as the critical mode ( 02 , wing forward bending) is 
concerned. However, sine feedback is unfavorable, in that the two modes which have the 
largest eigenvalue real part cannot be stabilized simultaneously. From figure 18, feedback 
of the cosine term of cyclic pitch is expected to work better than the sine term of cyclic 
pitch does. Consequently the feedback gain to the cosine term K^jcd is set to —0.03 for 
evaluation of this control system. 

Timewise calculation are made to see whether the feedback system just discussed is 
effective. Calculations are carried out without introducing rotor-blade or wing elasticity 
for the first five revolutions of the rotor rotation to get the air loading in an almost 
equilibrium state. The full equations of motion are employed after five revolutions; thus, 
the initial conditions for these equations are the deviations from the deflected equilibrium 
state. Figures 19 through 24 show the results for all the degrees of freedom both without 
and with feedback control. The abscissa is the number of rotor revolutions after elasticity 
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is introduced. The results for 100 revolutions, or approximately 15 sec, are given. With no 
feedback control, the wing-upward deflection appears to be converging (figure 19) because 
the lowest damped modes are dominated by generalized coordinates and as (wing- 
forward deflection and wing torsion). Wing-forward deflection (figure 20) and wing torsion 
(pylon-pitch angle; figure 21) show oscillation with almost constant amplitude with no 
feedback. Blade deformations do not appear to converge without control in spite of the 
fact that the lowest damped modes are not characterized by blade deflection. 

Once feedback control is introduced {w^ feedback with k^cd ~ -0.03) wing-forward 
and wing-torsional deflections are stabilized (Figs. 20 and 21). Stabilization is expected 
because these deflections are characterized by the two lowest damped eigenvalues and 
the real parts of each value are significantly reduced (stabilized) with the introduction of 
feedback control. The wing-upward deflection does not show an improvement because the 
eigenvalue, characterizing this motion (u; = 6.0 radjsec) is not affected by the introduction 
of feedback to the cosine term. The amplitude of blade-out-of-plane deflection becomes 
smaller since the eigenvalue for that mode (co = 62.3 radj sec) is changed by the control. 
The cyclic-pitch angle (figure 24) is sufficiently small so that there is additional control 
capability for stabilizing less stable flight conditions or for stabilizing oscillations of greater 
amplitude. From this time-history analysis, this control system has been proved to work as 
it should, and it would be expected to work as well under less stable operating conditions. 


CONCLUDING REMARKS 

Tilt rotor flutter control was investigated. The determination of the control law was 
based on a harmonic-balance method using a simple aerodynamic theory. It was confirmed 
that the stability of tilt-rotor motion is improved by using closed-loop, cyclic-pitch control. 
The validity of the control law determined in this way was demonstrated by a timewise 
calculation in which more sophisticated aerodynamics were considered. The selection of 
the feedback gain and the placement of the sensor were not based on any specific criterion. 
Establishing the optimal gain is left to future study. 
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APPENDIX A 


Expressions for , K*g, and 


Rotor-blade loading and are expressed in harmonic terms. 

Pen = Peo + P6s sin + p^c cos xl>n 

Pzn = PzO + Pzs sin V’n + Pzc COS 

Moreover the variables p#o» Peat Pec, Pzo, Pzs, and pzc are subdivided into pitch-angle 
and generalized-coordinate terms: 

2 3 

Pen =Peoo + Peoqu, 9jo + y~] Peoai^j 

j=\ i=\ 

2 2 

+ {pe SS 0s + E Pes^j.^js + ^^^PSsqj,,<ljc 
j=\ y=i 

3 3 

+ X] Peacij Oy + P^sa, ay) sin xjjn 
j=i y=i 

2 2 

{jPecc^c'^ 'y ^ Pecqjc Qjc “I” ^ ’'^ P0cq,,Qjs 


y=i y=i 

3 3 

+ ^ ^ Pecitj dy "h ^ ^ Pocaj gy) cos 
y=i y=i 

2 3 

Pzn =PzOO + PzOqjc QjO + PzOdj dy 

y=i y=i 

2 2 

■h {Pzss^s "I" ^ ^ Pzsqj, 9js “I" ^^Pzsqjr^jc 

y=i y=i 

3 3 

-h E Pzsdj^j H" E Pzsaj ) sin Xpfi 

j=l j=l 

2 2 

+ (P^cc^c + Pzcqj^Qjc + Pzcq^^ <]js 
j=l ; = 1 

3 3 

4" ^ ^ Pzcdj^j 4" ^ P2caj Xpyi 

j = l ;:=! 
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where 


1 - - oV^cC la [Bo - 4>) sin <f> - \pV^cCd cos 4> 
2v 1 - ^ 

1 ,r2 ^ 'S'WIJ . , 


2v/l -A/2 


pV^cCia sin (f> 


1 Tr2 /"t ^B2j . , 

s,n^ 

^ (>Vr cCio sin 

2v/r^ ^ vp 

- . =z pV^cCia sin (j) 

2yJ\-M^ ^ Vr 

i t /2 ^ ‘S'WSJ . , 

■ - / ~ -T7 ^PVRC<^la -T^ Sin <f> 
2v/l - M2 V^R 


pVpcCia sin<^ 


2n/T^" " " 

1 ,s2 /o ^BBj . , 

= r==^P^R^^lo‘~Tr^ 

2\/l - M2 "" V^j? 

—==^pV^cCia[Bo-<i>)coi4>~ -pV^cCds\n(l> 

2vl - M'' ^ 


2vT-A/2 


pV rCC la cos <f> 


2\/l - M2 


pv rCC la cos 4> 


2v/l - A/2 


pVRcCla cos 4> 
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2V1 - M2 


SB2j 
Vr 


pV^cC la COS <t> 


Pzsqjc - Vr 


Pzshj — 


Pzsaj — 


T - == pV^cCla cos <l> 

2 v/ 1 -M 2 ^ Vr 

1 .-2 ^ S’wsj , 

pVRcCia—r. COS(f> 


2V1 - M2 
1 




Pzcc ~ 


2v/l - M2 


pV^cCia COS (j) 


Pzcqjc — 


1 


2 n /1 - Af2 


Pzcqj, — 


Pzcdj — 


pVrcCic^^ COS (f> 

Vr 

^^.VjcC.^cos. 


2 vT-M 2 


Sw4j 

Vr 


pV rCC la cos (f> 


1 .,2 ^ Swsj , 

Substituting these relations into equation (29) yields 


ft — ^ + ^W2j^j) 


j=i 

2 


+ {^Blj^jO + ^B2j9j0 + f^BSj^js + ■^B 4 j 9 j 5 + ^BSj^jc + ^BQj^jc) 

y=i 


iV 

^Wlj ~ ~2 J [ ~ ^iPBshj + ^ZiPzOaj 

■!■ ^yi{~^P0cdj "I" ^Pzsij) 


~~ ^pi{^P6sdj "i" ^Pzcij) 


+ 2Nrirpeodj] dr 
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N 

^W2j ^ ~2 J [^t{~Pesaj + ^NyjPzOo) 

”H ^yi{ hpdcaj ^Pzsaj 2rJVpjP^Oo) 

“I" ^pi{~^^^yjP$00 ~~ ^POsUj hpzcaj)^ 

iv ' ^ 

^Blj ~ j (^^iPzOqjo ■(■ 2iVy,rp^Ogjo) 

N 

^B 2 j ~ ~2 J ’*' ^yi^P^sqj, ~ Npihpeaqj,^ dr 

N 

^BZj ~ ^ J ( ~ ^yihpocqjr ~ ^pi^Pzcqjc) dr 


Kh4, - 0 


K'b6j 


N 

~2 j f ~ ^iPOsqj, + ^yi{~^P0cqj. + ^Pzsqj,) 

“(■ Npii^jPzOO ~~ ^jPdOO ~~ ^Pzcqj, ~ f^P0sqj,)\ dr 

N 

~2 Jq I- ~ ^'P0«<ij-y 

"H ^yi{^jPzOO ^^jPOOO “I" ^Pzsqj,: hpffcqj,.) 

~~ ^pi [hpOsqjr. ^P^cqjc )] dr 

N 

Kg — ~~ I ( Kip$s 3 “I” Ky^rp^ss Kp{hpQas^ dr 

2 Jo 

N 

Kq — — I ( Kyihpffcc Kpirpj.cc) dr 

^ Jo 

Kq = N {Zipj.00 + Nrirpeoo) dr 

Jo 
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APPENDIX B 

Matrix Elements of Final Equations 


The following matrix elements for equations (112) through (114), with an 
for are given for: i = 1,2; j = 1,2; A: = 1,3; and / = 1,3: 



= 1 (t = 1 ~ 9) 

jp2 

^3t-2,/c4-6 

= Q\k 

^3i-l,k+6 

= Qik 

^Zi,k+6 

= Q\k 

jp2 

■^fc+6,3t-2 

= na{. 

jp2 

^k+6,Zi-l 

- 2 ^2t 

^k+6,3i 

- 2 ^3. 

^Zi-2,3j-2 

— — p 

- 

^3i-2,k+6 

— — P 

— ^Wlk 

Ipl 

— —P 

— ^B2j 

^3i-l,3i 

= -2Q 

^3i-l,k+6 

= Q3k “ ^\V2k 

jpl 

- 2Q 

^3t,3y 

— — p 

— ^B4j 

-^3t,fc+6 

— Q\k ~ ^W4k 

•^fc+6,3t-2 

= 

■^fc+6,3t-l 

= jAi, + NtlAt, - 

^k+6,3i 

= ^A\, - NUA\i - K%„ 

^k+6,l+6 

= lBNn{N,kNp, - NpiNpi,) 


G k T/-k 

W 3 l ~ ^Wll 

E”y(if not given above) 

Dtj(if not given above) 


jpO 

^3i-2,3i-2 

= 


jpO 

■^3t-l,3t-l 

= 

uf — 

jpO 

^3i-l,3j 

= 

_ ji 
^B3j 

17»0 

^3t-l,Jt+6 


_ ft 

^W3k 

jpO 

^3t,3j-l 

= 

— P 

^B5j 

J7»0 

^3i,3i 

= 

ujf - 

jpO 

^3t,A:+6 

= 

ft 

~^W5k 

rpO 

■^ik+6, 3i-2 


K^i 

jpO 

^it+6,3t-l 

= 

y(^*. + SI At, - 




17*0 

■^it+6,3i 

= 

- SI At, - 



l^k 

~ ^B6i 

jpO 

^>+6,/c+6 


i ~^W2k 


\ ~^W2k + ^wj 


Dzi-2,\ = Ibo 

Dzi-\,2 — Is 

Dzifi = Ic 

Dk+Q,l = Kq 

Dk+6,2 = Ks 

Dk+6,z = Kc 

= 0 

= 0 


exception 


n^ 4 .) 

n^A3^) 

j ¥= k 

j = k 
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APPENDIX C 

Vortex Lattice Method for Wing Aerodynamics 


Figure Cl shows the tilt-rotor-wing coordinate system for the bound- vortex calcula- 
tion. Only the right half wing is shown because the whole wing is assumed to be symmetric. 
The semispan wing is equally divided into k segments and the vortex S'li is trailed from 
the outer edge of the tth segment. Downwash at the tth control point v, is given as 

k 

Vi = ^ ^ 
j=\ 

where BijS'^j is the induced velocity on tth control point by the trailing vortex from the 
outer edge of jth segment. The wing’s mirror image and corresponding bound vortex 
connects the upper two trailing vortices. The coefficients Bij are given by 


+ Vbj ~ Vci ^ + ybj + Vci 

\/{yci - VbjV + Viya + VbjV + 

where Xc is the distance between the bound vortex and the control point and is equal to 
the semi-chord, y^j is the distance of jth segment outer edge from the center of whole wing, 
and yd is the distance of the ith control point from the center of the whole wing. Writing 
the boundary condition (i.e., downwash on the tth control point) as 6i, the equation which 
connects the intensity of each trailing vortex and the boundary conditions is given by 


47r(y6j - yd) 


{B) 


fS'n \ 




\bkJ 


(B) = 


/Bn 

Bn 

... Bik 

B21 

B22 

. • • B^k 

VBjki 

Bk 2 

• • • Bkk 


The intensity of the trailing vortex is given by 


\ 

H2 


. . . 



UJ 
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The intensity of the bound vortex in the jth section 'yty is 


Ibk = Hk 

lbk-i = lk+^lk-x t = 1,2,. . . A: - 1 

Finally, the lift distribution in tth section /, is 

U = pVibi 

To determine the direction of the lift, the induced velocity on the bound vortex v\,i needs to 
be determined. The lift is caused by the trailing vortex only, since a straight and unswept 
wing is considered. The equation for the induced velocity is 

47T " V Vbi - Vci Vbj + Vex J 
The backward inclination of lift (j>w is given by 

(f>xjii —VbifV 

^ -A / S'lj \ 

\ Vbj - ya Vbj + Vex ) 
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Figure Cl.- Vortex lattice geometry for wing aerodynamics. 
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TABLE 1.- SPECIFICATION OF BASELINE TILT ROTOR MODEL. 


ROTOR 

number of blades, N 

3 

radius, R 

3.81 m 

cutout, Rc 

0.15R 

chord , c 

0.356 m 

solidity, a 

0.076 

lift slope, Cia 

5.73 

drag coefficient, Cj, 

0.0125 

rotational angular velocity, D 

40.4 rad/sed (386.0 rpm) 

blade-flapping inertia. Is 

176.1 kg m2 

blade mass, Mb 

42.9 kg 

WING 

semispan, L 

5.08 m 

chord, 

1.58 m 

geometric angle of attack, aw 

2.7° 

lift slope, Cia 

5.73 

drag coefficient, cj, 

0.004 

pitching moment coefficient. Cm 

-0.005 

position of elastic axis from leading edge, Cg 

0.26 

position of center of gravity from leading edge, Cg 

0.26 

PYLON 

mass, Mp 

644.0 kg 

mast height, h 

1.3 m 

pitching mass moment of inertia, Ipp 

257.6 kg • m2 

rolling mass moment of inertia, Ipr 

57.5 kg • m2 

yawing mass moment of inertia, Ipy 

223.4 kg • m2 
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Figure 1.- Tilt-rotor model. 


X 



Figure 2.- Coordinate system and pylon rotations. 
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a.- Wing. 



ROTATIONAL 

PLANE 



b.- Blade. 


Figure 3.- Wing and blade deformation. 
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Figure 4.- Rotor azimuth angle and pylon rotation. 
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Figure 7.- Blade structural properties. 
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a.- Mass and moment of inertia. 
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b.- Stiffness. 

Figure 8.- Baseline-wing structural properties. 
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Figure 9.- Pitch angle and thrust necessary for cruise flight. 
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Figure 10.- Real part of eigenvalues for baseline-wing configuration. 
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Figure 11.- First blade mode shape; u>i = 37.6 rad /sec. 



Figure 12.- Second blade mode shape; W 2 







Figure 13.- First wing mode shape; Wuii = 6.19 rad/sec. 



Figure 14.- Second wing mode shape; Wu;2 = 11.2 rad/sec. 
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Figure 15.- Third wing mode shape; Wu ,3 = 18.1 rad/sec. 



Figure 16.- Real part of eigenvalues for one eighth stiffness wing configuration. 
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Figure 17.- Eigenvectors without feedback for lowest damp>ed modes; V (QR = 1.5. 
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FEEDBACK GAIN 

a.- Cosine cyclic feedback; K^^d- 


Figure 18.- Real part of eigenvalues with wing feedback; VfflR = 1.5. 
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WING UPWARD 
DEFLECTION, m 


0 10 20 30 40 50 60 70 80 90 100 

REVOLUTION 


a.- No feedback control. 



REVOLUTION 


b.- Feedback control; — —0.03. 
Figure 19.- Wingtip upward deflection. 
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REVOLUTION 


a.- No feedback control. 



REVOLUTION 


b.- Feedback control; Kujcd = —0.03. 
Figure 20.- Wingtip forward deflection. 
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PYLON PITCH PYLON PITCH 

ANGLE, radian ANGLE, radian 



REVOLUTION 

a.- No feedback control. 



REVOLUTION 

b.- Feedback control; Kwcd = —0.03. 
Figure 21.- Pylon-pitch angle. 
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BLADE IN-PLANE BLADE IN-PLANE 

DEFLECTION, m DEFLECTION, m 


.2 



•2 I I I I 1 I I I I I 1 

0 10 20 30 40 50 60 70 80 90 100 

REVOLUTION 


a.- No feedback control. 


•2 r 



0 10 20 30 40 50 60 70 80 90 100 

REVOLUTION 


b.- Feedback control; K^jcd = —0.03. 


Figure 22.- First-blade inplane deflection. 
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Figure 24.- Cosine cyclic-control pitch used for feedback; Kwcd = -0.03. 
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